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1 Introduction 



The Hanbury-Brown-Twiss (HBT) effect [0 originated in astronomy where one uses the 
interference pattern of the photons to learn about the size of the photon emission region, i.e. 
the size of the particular star, which is emitting the light. The effect can be described as an 
enhancement of the two-particle correlation function that occur when the two particles are 
identical bosons and have very similar energy-momenta. A well-known formula to relate 
the two-particle correlation function (in four-momenta pj,j = 1,2 with q = pi — P2) to the 
space-time density distribution, p, of (chaotic) emission sources is, 

+ (1) 



daida2 

where TZ is the normalised Fourier transform of the source density 

J p{x)dx expjiqx) 

This quantity is often, without very convincing reasons, parametrised in terms of a "source 
radius" R and a "chaoticity parameter" A, 

n{q) = Xexp{-R^Q^) (3) 

with = — g^. The source radii obtained by this parametrisation tend to be similar in 
all hadronic interations (we exclude heavy ion interactions where the extensions of nuclear 
targets and probes will influence the result), with R ~ 0.5 — 1 fm, but the chaoticity varies 
rather much depending upon the particular data sample and the method of the fit. At 
present the knowledge of higher-order correlations is still limited in the experimental data, 
although in principle there should be such correlations. 

The HBT effect, being of an intrinsically quantum mechanical nature, is not easily 
included in the event generator programs used in high energy physics. Such simulation 
programs, like HERWIG (based upon the Webber-Marchesini parton cascades and ending 
by cluster fragmentation) and JETSET ||^ (based upon the Lund Model string dynamics [0) 
are built in accordance with classical stochastical processes, i.e. they produce a probability 
weight for an event without any quantum mechanical interference effects. 

Sjostrand has introduced a clever device as a subroutine to JETSET, in which the HBT 
effect is simulated as a mean field potential attraction between identical bosons [Q. Thus, 
given a set of energy-momentum vectors of identical bosons, pi, ■ ■ ■ ,Pn, generated without 
any HBT effect, it is possible to reshuffle the set into another set where each pair on the 
average has been moved relatively closer to show a (chosen) HBT-correlation, while still 
keeping to energy- momentum conservation for the whole event. 

In this paper we will develop a method deviced in to provide the Lund Model with a 
quantum mechanical interpretation. In particular there will be a production matrix element 
with well-defined phases. This will then be used to make a model of the HBT effect. Although 
this model stems from different considerations it will nevertheless contain predictions which 
are similar to those in the ordinary approach giving Eq (Q). 
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In the first section we survey those features of the Lund model, that are necessary for 
the following. In the second section we will exhibit the general n-particle HBT-correlations 
in the model. The resulting expressions contain a sum of in general n! terms, i.e. it is 
of exponential type from a computational point of view. It is possible to subdivide the 
expressions in accordance with the group structure of the permutation group. Although 
the higher order terms provide small contributions in general the computing times are still 
forbidding. In order to speed up the calculations we introduce instead the notion of links 
between the particles and it is in this way possible to obtain expressions of power type from 
a calculational point of view, which are perfectly tractable in a computer. In the last section 
we exhibit a set of results both in order to show the workings of the model and to provide 
predictions for experiments. 

We will in this paper be satisfied to treat only two-jet events, i.e. we will neglect hard 
gluon radiation and we will come back to HBT effects in gluon events in another publication. 



2 Some Properties of the Lund Model 

Within the framework of perturbative QCD it is possible to obtain many useful formulas 
but all the results are expressed in a partonic language. In order to be able to compare to 
the hadronic distributions, which are observed in the experimental setups, it is necessary to 
supplement the perturbative results with a fragmentation process. We will in this paper be 
concerned with the Lund string model and we start with a brief introduction to its main 
properties. 

In the string model the confining colour field is approximated by a massless relativistic 
string. The endpoints of the string are identified with quark and antiquark properties while 
the gluons are assumed to behave as transverse excitations on the string. The string can 
break up into smaller pieces by the production of gg-pairs (i.e. new endpoints). Such a pair 
will immediately start to separate because of the string tension, which in the rest frame of a 
string segment corresponds to a constant force k; phenomenologically k ^ 1 GeV/ fm. Final 
state mesons are formed from a q and a q from adjacent vertices, as shown in Fig |l|. 

Each breakup vertex will separate the string into two causally disconnected parts. From 
the causality, together with Lorentz covariance and straightforward kinematics, it is possible 
to derive a unique breakup rule for the string by means of (semi) classical arguments 

The unique breakup rule results in the following probability for a string to decay into 
hadrons (pi, . . . ,p„). 



dP{pi, ...,pr, 



6iY.Pj-Ptot)expi-bA) (4) 



where A is the area of the breakup region as indicated in Fig ^ and N and b are two 
parameters. 

The similarity of the result to Fermi's Golden Rule for the probability of a quantum 
mechanical transition, i.e. the size of the final state phase space multiplied by the square of 
a matrix element \Ai\'^ expressed by the exponential area suppression, provides a reasonable 
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Figure 1: The decay of a Lund Model string. 



starting point to try to derive a corresponding quantum mechanical process. There is at 
least two possible mechanisms, viz. a quantum mechanical tunneling process a la Schwinger 
and/or the possible relationship to the Wilson loop operators in a gauge field theory. We 
will find that they provide very similar answers to the problem. 

Starting with the tunneling arguments, we note that while a massless g'^-pair without 
transverse momentum can be produced in a pointlike way anywhere along the string, a 
massive pair or a pair with transverse momentum must classically be produced at a distance 
so that the string energy between them can be used to fulfil energy- momentum conservation. 
If the transverse momentum is conserved in the production process, i.e. the qq with masses ii 
obtain ±^x, respectively, then the pair may classically be realised at a distance 5x — 2^j^/k, 



The probability for a quantum mechanical fluctuation of a pair, occurring with /x^ at 
the (space-like) distance 5x, is in a force-free region given by the free Feynman propagator 
squared: 



A corresponding quantum mechanical tunneling process in a constant force fleld will accord- 
ing to WKB methods give 



The difference is that in the force-free case we obtain an exponential suppression 4//^//« 
but when the constant force pulls the pair apart we obtain the somewhat smaller suppres- 
sion TTyU^/zt. Besides the mass suppression (which phenomenologically will suppress strange 
quark-pairs with a factor of ~ 0.3 compared to "massless" up and down flavored pairs) we 








3 



obtain the transverse momentum Gaussian suppression 



exp(-^kl) with 2a^ = - (7) 

The value of a as used in JETSET is a bit larger than the result in Eq but this can be 
understood as an effect of soft gluon generation along the string. The transverse momentum 
of a hadron produced in the Lund Model is then the sum of the transverse momenta of its 
constituents. 

We may use the elementary result in Eq to calculate the persistence probability of 
the vacuum, V, as it is defined in It is the probability that the no-particle vacuum will 
not break up, owing to pair-production, during the time T over a transverse region A^, when 
a constant force k is apphed along the longitudinal x-direction over a region L: 



V= n (l-P(/i^))=exp 

t£{0,T),x€{0,L),k±,s,f 



E log(i-p) 

\t,x,k±,s,f 



We have then assumed that the field couples to (fermion) pairs with spin s and flavors 
/ and we sum over all possibilities for the production. As each pair needs a longitudinal 
size Sx = 2fi±/K, and, according to Heisenberg's indeterminacy relation, will live during a 
timespan 27r/2/i^ there is at most k,LT/2tt pairs possible over the space-time region LT. The 
transverse momentum summation can be done by Gaussian integrals from an expansion of 
log(l — P) and the introduction of the well-known number of waves available in a transverse 
region A±: {A±/{2TTy)d'^k±. In this way we obtain for the persistence probability 



°° 1 / , ,2 ' 



V = expi-^LTA^U) with n = ^V^exp (9) 



=1 



where nf,ns is the number of flavor and spin states. 

There are two remarks to this result. Firstly, although the method to treat the inte- 
gration over time and longitudinal space, by close-packing reasonably sized boxes, may not 
seem convincing the final formula coincides with the one obtained by Schwinger [^, for 
the case of a constant electric field S. Then n is identified with the force of the charges in 
the external field, i.e. n — > eS. 

Secondly, the result is in evident agreement with the formula for the decay of the Lund 
string in Eq if we identify LT with the (coordinate space) area size A. In this way we 
also obtain the result that the parameter b is 

b = k'^A^U (10) 

i.e. it corresponds to the transverse size of the (constant) force field, which we have modelled 
by the string. The quantity 11 is l/(127r) for two massless spin 1/2-fiavors. 

The second quantum mechanical approach is to note that a final state hadron stems 
from a q from one vertex j and a q from the adjoining vertex j + 1. In order to keep to 
gauge invariance it is then necessary that the production matrix element contains at least a 
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gauge connector between the vertices: exp{i J^~^^ gA^dx^), where g is the charge and the 
gauge field. Consequently the total production matrix element must contain a Wilson loop 
operator: 

M = exp(z <j) gA^dx^,) (11) 

with the integration around the region A (note that the field is singular along the border line 
and we are therefore not allowed to distort the integration contour inwards). The operator 
in Eq ([TTl) was predicted (and inside lattice gauge calculations also found) to behave as 

M = ex.p{i^A) (12) 

with the real part of ^, Re{^) = k. In the present situation where the force field region 
decays we expect an imaginary part, corresponding to the pair production rate according to 
the well-known Kramers-Kronig |]lT| relationship for the dielectricity in matter, in this case 



the QCD vacuum. 

The two interpretations of the area law, i.e. the Schwinger tunneling in Eq (y) and the 
Wilson loop operator result in Eq (|l^) can be related if we note that according to Gauss' 
Law the integral over the extension of the force field should correspond to the charge. For 
a thin string we should then obtain for the area falloff rate b oc k^A± oc na. Although 
Gauss' law is more complicated for a non-abelian field with triplet and octet color-charges 
and similarly octet fields it is possible to make a case for an identification of the parameter 
6 as 

b = ^ (13) 
12 ^ ^ 

which is what we should expect from the expected imaginary part of the dielectricity in 
Eq ([121) . a = 3g'^/{4:7i) is then the effective QCD coupling, including the color factors. 
The result is also phenomenologically supported if we consider a partonic cascade down 
to a certain transverse momentum cutoff k±c and then use the Lund model hadronisation 
formulas to obtain the observed properties of the final state. In that way we may determine 
the parameters in the model as functions of the partonic cascade cutoff. A remarkably good 
fit to the 6-parameter is given by C/ loglk^^/A"^) with C given by Eq (p!3D and A ^ .5 GeV 
H^, according to the QCD coupling. 

Independently of the precise identification of b, we obtain a possible matrix element 
from Eq (0) 

M = exp{iK - b/2)A (14) 

which not only will provide us with the Lund decay probability in Eq (^, but also can be 
used in accordance with to provide a model for the Hanbury-Brown-Twiss effect (often 
called the Bose-Einstein effect) for the correlations among identical bosons. 



3 A Model for Bose-Einstein Correlations 

We will from now on work in energy-momentum space in agreement with the usual treatment 
of the Lund model formulas. Further we will make use of a lightcone metric with p± = e±pe 
where i denotes the longitudinal direction along the string. The two metrics differ by a 



5 



factor of two, i.e. 2dedpi = dp^dp^. Note in particular that compared to the considerations 
in the earher section this means that the area A — * and the parameter b — > 6/(2k^). 

We now consider a final state containing (among possibly a lot of other stuff) n identical 
bosons. There are n\ ways to produce such a state, each corresponding to a different per- 
mutation of the particles. According to quantum mechanics the transition matrix element is 
to be symmetrised with respect to exchange of identical bosons. This leads to the following 
general expression for the production amplitude 

M = Y.Mv (15) 

V 

where the sum goes over all possible permutations V of the identical bosons. The cross 
section will then contain the square of the symmetrised amplitude Ai 

lA^I-^BlA^H-d+E i^T/l^-' )) (16) 

•p -pi^-p \-Mv\ + I-M-P'l 

The phenomenological models used to describe the hadronisation process are formulated 
in a probalistic language (and not in an amplitude based language). This implies that 
interference between different ways to produce identical bosons is not included. In this case 
the probability for producing the state is 

\M\' = Y.\Mv? (17) 

T 

instead of the probability in Eq (|^). Comparing Eq (0) and Eq (|1^ it is seen that a 
particular production configuration leading to the final state V can be produced according 
to a probalistic scheme and that the quantum mechanical interference from production of 
identical bosons can be incorporated by weighting the produced event with 

The outer sum in Eq (|l^) is as usual taken care of by generating many events. 

In order to see the main feature of symmetrising the hadron production amplitude 
in the Lund Model we consider Fig H, in which two of the produced hadrons, denoted 
(1,2), are assumed to be identical bosons and the state in between them is denoted I. We 
note that there are two different ways to produce the entire state corresponding to the 
production configurations (. . . , 1, /, 2, . . .) and (. . . , 2, /, 1, . . .), i.e. to exchanging the two 
identical bosons. The two production configurations are shown in the figure and the main 
observation is that they in general correspond to different areas! 

The area difference, AA, depends only on the energy momentum vectors pi-,P2 and p/, 
but can in a dimensionless and intuitively useful way be written 

= 5p5x (19) 

where 5p = p2—pi and Sx = {St; 0, 0, 6z) is a reasonable estimate of the space-time difference, 
along the surface area, between the production points of the two identical bosons. We note 
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Figure 2: The two possible ways to produce the entire state when two of the bosons are 
identical. 



that the space-time difference Sx is always space-like. In Fig ^ Sx, for the two production 
configurations, is indicated by arrows, together with open circles showing the corresponding 
production points. 

We go on to consider the effects of transverse momentum generation in the gg-vertices. 
First we note that the total transverse momenta of the substate 1, /, 2 in Fig ^ stem from the 
q and q generated at the two surrounding vertices, A and B. This is, owing to momentum 
conservation, fixed by the properties of the hadrons generated outside of the substate. Using 
this we find that there is a unique way to change the transverse momenta in the vertices sur- 
rounding the intermediate state / such that every hadron has the same transverse momenta 
in both production configurations. 

Suppose as an example that we have generated ik^y^ in the vertex A and ik^^ in the 
vertex B (i.e. so that — k_|_yi and k_|_B defines the substate). Then to conserve the transverse 
momenta of the observed hadrons when changing production configuration from (1, /, 2) to 
(2, /, 1) it is necessary to change the generation of transverse momenta in the two vertices 
surrounding I as follows (in an easily understood notation): 

± k^j ^ ±{k^A + kxB - Ki) (20) 
±k'^^ ^ ±{k±A + kxB - 

This means that exchanging two bosons with different transverse momenta will result in a 
change in the amplitude as given by Eq (0) for some of the vertices. 

From the amplitudes in Eq (|T^ and Eq (^ we get that the weight in the Lund Model 
can be written 

w = i+y] — 7 , (21) 

-^-cosh ^^ + ^^' 
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where A denotes the difference with respect to the configurations V and V and the sum of 
p]_^ is over all vertices. We have introduced cTp^ as the width of the transverse momenta for 
the generated hadrons, (i.e. = 2(T^). 

Using Eq (|l^) for a single pair exchange one sees that the area difference is, for small 
governed by the distance between the production points and that AA increases quickly with 
this distance. We also note that AA vanishes with the four-momentum difference and that 
the contribution to the weight from a given configuration, V , vanishes fast with increasing 
area difference AA. From these considerations it is obvious that only exchanges of pairs with 
a small 5p and a small 5x will give a contribution to the weight. In this way it is possible to 
relate to the ordinary way to interpret the HBT effect, cf. Eq (13). 



It is straightforward to generalise Eq (|I9|) to higher order correlations. One notes in 
particular that the area difference does not vanish if more than two identical bosons are 
permuted and only two of the bosons have identical four-momenta. 



Models for BE correlations has been suggested, e.g. |T3|, with similar weight functions, 
but it is important to note that the weight in our model has a scale both for the argument to 
the cos-function as well as for the function which works as a cut-off for large 5p and 5x (in our 
case a cosh- function). Further the two scales in our model are different and well-defined, at 
least phenomenologically. We will come back to the influence of the two scales in Section ^ 



4 MC Implementation 

To calculate the weight for a general event, with multiplicity n, one has to go through 
ni!n2!...n7v! — 1 possible production configurations, where is the number of particles of 
type i and there are N different kinds of bosons. For a general e"'"e~-event at 90 GeV this is 
not possible from a computational point of view. 

We know however that the vast majority of configurations will give large area differences 
and they will therefore not contribute to the weight. One of our aims with this work has 
been to find a way to approximate the sum in Eq ( PT| ) with a sum over configurations with 
significant contributions to the weight. From basic group theory we know that every group 
can be partitioned into its classes. Let 11111 ... 1 denote the class containing only the identity 
element, where all particles are unchanged, mim2lll ... 1 denote the class of group-elements 
where mi particles are cyclically permuted, m2 other particles are cyclically permuted and 
the rest are unchanged, and so on. We can define the order, /c, of a class as k = Y^iij^i — !)• 
The useful feature of this ordering of classes is that for all group-elements contained in order 
k the minimum of the summed size (in positions) of the cyclically permuted clusters , Ar, is 
fc, i.e. the minimum length over which particles are moved increases with the order. From 
the discussion at the end of Section ^ it is then obvious that the contribution to the weight 
from a configuration will decrease with its order. All classes up to order 4 are shown in 
Table |. 

We have found that for essentially all events the weight does not change when including 
the fifth order. But we have also found that lots of lower-order configurations give no 
contributions to the weight. This is not acceptable when taking computing time into account 
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Order 





1 


2 


3 


4 


Classes 


11111.. .1 


21111... 1 


22111. ..1 
31111. ..1 


22211. ..1 
32111. ..1 
41111. ..1 


22221... 1 
32211... 1 
33111... 1 
51111... 1 
42111... 1 


Ar 





> 1 


> 2 


> 3 


> 4 



Table 1: The classes of the permutation group order by order up to the fourth order. Ar is 
the minimum length over which particles are permuted. 

and we have therefore abandoned using a cut in order. 

In this work we have instead approximated the sum in Eq (0) with a sum over con- 
figurations of all orders with significant contributions to the weight. This has been done by 
introducing exchange-links between particles. We have only taken into account interference 
with configurations where all particles are produced in positions from which there is a link 
to a particle's original production position. Defining an exchange matrix, C, as follows 

^ f 1 if there is a link between particles i and j. 
1 otherwise 

one gets a simple representation of the configurations to be considered. If all elements in C 
are 1, it corresponds to considering all n! permutations, while only the original configuration 
is considered if C is the identity matrix. 

AA for a pair exchange increases, as previously discussed, with the four- vector difference 
6p and with the size of the state inbetween. Since we know from Eq (^) that the contribution 
to the weight for a given configuration vanishes fast with increasing area difference A A, it 
is useful to introduce the concept of linksize, defined below as the invariant four-momentum 
difference together with the invariant mass of the particles produced in between the pair 
(in rank). By only accepting links between particles if the size of the link between them is 
smaller than some cut-off linksize, Sc, we get a prescription for the exchange matrix of an 
event. In this way, by specifying the allowed two-particle exchanges, we get, to all orders, 
which configurations to take into account. We have found that for a given Sc one includes all 
configurations that provide a contribution larger than some e to the weight. Taken together 
this means that we get all the important contributions to the weight if we chose 6c so large 
that the neglected terms smaller than e give a negligable change for every weight. 

We have used a cut-off linksize such that there is a link between identical bosons if one 
of the following conditions is fulfilled. 

• = -{p^ - p,r < Qla. ^ ^ GeV\ 

• the invariant mass of the particles produced inbetween (along the string) the pair is 
less than m^„^ ~ (20 GeVf. 
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Including links larger than this give no contribution to the weight for essentially all events. 
There are a few special events for which the weights have not converged with this 5^- They are 
very rare and have in common that they have a cluster of particles such that exchanging any 
pair in the cluster will give a large area-difference, but there are cyclic permutations which 
give a small area-difference. Increasing 5c to include these configurations give no noticeable 
effect in any observable known to us (except the computing time in the simulation!). 



Including decays 

A large fraction of all final-state bosons stem from decays of short-lived resonances with 
lifetimes comparable to the time scale in string decay. Therefore they may contribute to 
the Bose-Einstein effect. To include their decay amplitudes and phase space factors and 
symmetrise the total amplitude is very difficult and it is furthermore not known how to do 
that in a model-consistent way. We have included resonance decays in the following simple 
way 



Particles with width larger than Tmin are assumed to decay before Bose-Einstein 
symmetrisation sets in and the matrix elements are evaluated with their decay 
products regarded as being produced directly, ordered in rank. We have used 



The signal in the two-particle correlation function goes down very much if we neglect all 
the pions from resonance decays when symmetrising the amplitudes. But our signals are 
fairly independent of T„dn as long as it is small enough for the p's to decay before the 
symmetrisation. 

An elaborate discussion on the treatment of resonances in connection with BE correla- 



tions can be found in 14 



5 Results 

In our simulations we have used the Lund string model 0] implemented in the JETSET MC 
^ to hadronize qq-pairs (i.e. no gluons are considered). The MC implementation of our 
model is available from the authors. 



5.1 The Weight Distribution and Two-particle Correlations 

The majority of the weights are close to and centered around unity, as seen in Fig |[ There is 
however a tail of weights far away from unity in both directions. The tail of positive weights 
is shown as an insert and the distribution looks like a power. However if we subdivide the 
events into sets with similar number of links and study the weight distributions for these sets 
separately, we find that the weight distribution for each set is basically Gaussian. The width 
of these Gaussians increase with the number of links in the corresponding set, as shown in 
Fig ^. The power like behaviour of the weight distribution is therefore merely a consequence 
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Figure 3: The distribution of Bose-Einstein weights for two-jet states in JETSET. The tail 
of positive weights is shown in the insert. 



-2 




10 100 

¥eight(w) 

Figure 4: The distribution of Bose-Einstein weights for two-jet events subdivided into sets 
with different number of links, rti. Two samples for 3 < n/ < 5 and 10 < ri; < 12 are plotted. 

of summing over events with different number of links. It should be emphasized that the 
negative weights only are a technical problem. Summing over many events results in positive 
probabilites for all physical observables, which is obvious from Eq (|T^). 

We have taken the ratio of the two-particle probability density of pions, p2, with and 
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without BE weights apphed as the two-particle correlation function, R2, i.e. 



^2(^1,^2) 



p2w{Pl,P2) 
P2{PI,P2) 



(22) 



where the w denotes weighted distributions. 



As discussed in connection with Eq (^9]) the correlation length in Q depends inversely 
on the (space-like) distance between the production points of the identical bosons and the 
Bose-Einstein correlation length, that is dynamically implemented, in this model can most 
easily be described as the flavour compensation length, i.e. the region over which a particular 
flavour is neutralised. Identically charged particles cannot be produced as neighbours along 
the string in the Lund model while neutral particles can. This implies that identically 
charged pions which always must have a non-vanishing state inbetween will have a more 
narrow correlation distribution in Q compared to neutral pions. This has been found as 
can be seen in Fig |^ where the correlation distributions for pairs of particles used in the 
symmetrisation are shown. The correlation functions have been normalised to unity in the 
region 1.0 < Q < 2.5. The correlation distribution for charged pions can be approximated 
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Figure 5: The ratio R{Q) of the number of pairs with invariant relative four-momentum Q 
with and without Bose-Einstein weights applied. The sample consists of the initial particles 
used in the symmetrisation. 

by the LUBOEI algorithm [^] with radii ~ 1 fm and A ~ 0.8 as input parameters (Note 
that the input parameters are not exactly reproduced in the resulting correlation function). 
This is in reasonable agreement with the LEP experiments, which measure sizes of the 
order of 0.5 — 1 fm, |l^, |16|, [l^, even though they tend towards values smaller than 1 fm. 
The correlation function depends on the daughters of resonances and especially the decay 
products of rj' play a large role. The production rate of rj' used in JETSET was questioned in 
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, in connection with BE correlations. We have used reduced production rates for rj and t]' 



by setting the extra supression factors in JETSET to 0.7 and 0.2 respectively, in accordance 



with the DELPHI tuning |T8[. For a more elaborate quantative comparison with data our 
BE Monte Carlo has to be tuned further and the resulting events has to be subjected to the 
same corrections as in the experimental analysis. 

The general findings for the parameter dependence of the weight function in Eq (^Tj) is 
that due to the smallness of b as compared to I/k, it is for most of the terms in the sum 
the decrease of the cos-function with increasing AA that governs the behaviour. For larger 

it is the transverse momentum contributions to the cosh-function which takes over to 
damp the contribution to the weight. Note that the argument of the cos-function contains 
K as the basic scale and that the transverse momentum contributions also are governed 
essentially by k, (see Eq (|^)). Going over to correlation functions we find as expected from 
the conclusions for the weight function that the correlation function is not affected when the 
6-parameter is changed ±20%. The slope of the correlation function for small Q values and 
therefore the correlation length is very sensitive to k and it is also sensitive to the width of 
the transverse momenta. The transverse momentum generation acts as noise in the model so 
that all weights approach unity and consequently all correlations vanish with increasing o"p^. 
It is however this noise which makes the weight calculations tractable. Consequently, the 
main parameter is the string tension, k, in this model for Bose-Einstein correlation weights 
as well as for the correlation length. 

Since the BE weights are depending on the space-like distance between the production 
points we have studied the two-particle correlation function as a function of the invariant 
space-time distance Ax = ^ —bx^ where bx = {5t; 0, 0, 6z) as defined in Eq (|19[) . In Fig ^ 
R2{Ax), which has been normalised to unity in the region 1.0 < Q < 2.5, is plotted. The 
figure illustrates that the effect of the Bose-Einstein symmetrisation , i.e. to pack identical 
bosons closer together in phase-space, is manifest up to production point separations of about 
0.7 fm. It should however be noted that many configurations where pairs are exchanged 
over significantly larger distances give significant contributions to the weight. 

We have also found that the higher order contributions to the sum in Eq (^) is of 
importance for the two-particle correlations. That is using more than two-particle exchanges 
when calculating the weights does not only affect the weight distribution but also the two- 
particle correlation function, R2{Q)- 

In heavy-ion collision experiments one has found that the extracted correlation length 
has an approximate l/y/rn^ dependence [^]. This is in agreement with hydrodynamical 
models describing the source evolution in heavy-ion collisions. Recently a similar rrit de- 



pendence has been found for hadronic decays in e'^e annihilation at LEP |2^. In the 



Lund Model the average space-like distance between pairs of identical pions increases with 
rrit and one would therefore not expect a correlation length which falls off with rrit. For 
initially produced particles we get a correlation length which is essentially independent of 
rrit. However when analysing all final particles we find for increasing rrit that the correlation 



length falls off and that the A parameter increases, as in [^. From this we conclude that 
the observed rut dependence of the correlation length in data is, in our model, compatible 
with the vanishing of contributions from decay products with increasing m^. 
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Figure 6: The ratio R2{Ax) of the number of pairs with invariant relative distance Ax with 
and without Bose-Einstein weights applied. The sample consists of the initial particles used 
in the symmetrisation. 



5.2 Residual Bose-Einstein correlations 

Bose-Einstein correlations acting between identical bosons may have significant indirect ef- 
fects on the phase space for pairs of non-identical bosons. We have studied mass distributions 
of TT+TT^ systems to see how our model affects systems of unlike charged pions. Many analyses 
use TT'^7c~ distributions to quantify the Bose-Einstein correlations, using the unlike-charged 
distributions as reference samples with which to compare the like-charged pion distributions. 
We have found that the assumption that the two-particle phase space densities for vr+vr" sys- 
tems are relatively unaffected by Bose-Einstein symmetrisation is fairly good. Taking the 
ratio of the vr+vr" mass distributions with and without Bose-Einstein symmetrisation applied 
gives that the mass distribution is not altered much by the symmetrisation, and that the 
effect is smaller than 5% in the entire mass range. 

It has however been observed experimentally that the Breit-Wigner shape for oppositely 



charged pions from the decay of the p resonance [^, |21|, ^ is distorted. We have therefore 
analysed tt'^tt~ distributions when the pair comes from the decay of a These vr+vr" mass 
distributions, with and without BE weights applied, are shown together with the difference 
of the two in Fig ^ From the difference it is clearly seen that the weighting depletes the 
region around the p mass and shifts the masses towards lower values as well as it slightly 
increases the width of the distribution. The figure clearly shows the potential of our model 
to affect the mass spectrum of the p°. 
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Figure 7: p° meson mass shift induced by the Bose-Einstein correlations in our model. The 
solid curve shows the Breit-Wigner as generated by JETSET while the dashed curve is ob- 
tained after applying BE weights to the events. The curve with error bars is the difference of 
the two (dashed- solid). The areas under the two mass distributions are normalised to unity 
in the shown mass range. 



5.3 Three-particle correlations 

The existence of higher order dynamical correlations, which are not a consequence of two- 
particle correlations, is of importance for the understanding of BE correlations. There are 
very few experimental studies of genuine three-particle correlations, mainly because of the 
problem of subtracting the consequenses of two-particle correlations and the need for high 
statistics of large multiplicity events. Genuine short-range three-particle correlations have 
been observed in e"'"e~ annihilations by the DELPHI experiment. They conclude that they 



can be explained as a higher order Bose-Einstein effect |23 



To reduce problems with pseudo-correlations due to the summation of events with dif- 
ferent multiplicities we have used three-particle densities normalised to unity separately for 
every multiplicity in the following way 



~(a,6,c)/ N 



~{na,nf,,na) 



(Pl,P2,P3) 



n>8 



~('na,ni,,nc) / \ 
P3 {Pl,P2,P3] 



nairib - Sab){nc - Sac - Sbc) 0-{na,ni,nc) dpidp2dp3 



(23) 



(24) 



where n is the charged multiplicity, cr„^,nj,nc is the semi-inclusive cross section for events with 
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rii particles of species i, and 



PK, n,, = "^^"-"-■"-^ (25) 

Z^na,ni,,nc (ua ,ni, ,nc) 

We have aimed to study the genuine normahsed three particle correlation function, R^, 
defined as 

^3 = [P3(Pl,P2,P3) - p2{Pl,P2)pl{P3) " p2iPl, Ps) Pl{P2) " p2{P2, Ps) Pl{Pl) 

+2pi{pi)pi{p2)pl{p3)]/{pl{Pl)pl{P2)pl{P3)) + 1 (26) 



where we have used an abbreviated notation for the pa from Eq (|23D, and pi and p2 are the 
corresponding one- and two-particle densities, normalised in accordance with Eq ( P5[ ) and 
Eq (0). R3 is equal to one if all three-particle correlations are consequenses of two-particle 
correlations. 

In order to calculate the P2P1 and pipipi terms in Eq (^) the common experimental 
procedure is to mix tracks from different events. Using a mixing procedure in our model 
means weighting triplets of particles with products of event weights. This results in large 
statistical fluctuations and to get them under control, with our event weights, requires gener- 
ation of very many events. We have therefore taken another approach, in order to minimise 
the computing time. We have used combinations of charged pions in the following way to 
approximate Eq (PB| ) 

5 _ PSW "^[PSW PS ) /r.'TN 

^3 = -(±,±,±) (^') 

P3 

where w, as previously, denotes weighted distributions. There are a couple of things to 
note in connection with Eq (pTl). If there are genuine positive three-particle correlations for 



(+ H — ) and ( h) combinations, as observed by the DELPHI collaboration |^ they will 

if they come from BE symmetrisation contribute to the R3 in Eq (|27|), but they will reduce 
the signal. Secondly, we note that there is a possible bias from two-particle correlations from 
(H — ) combinations but that it is small as discussed previously. We also note that using the 
normalisation in Eq (p^ reduces problems with contributions from like- and unlike-charge 
combinations having different multiplicity dependence. It should also be observed that the 
i?3 in Eq (|^) can be studied experimentally since getting the ps^j's of course is achieved by 
analysing single events and the p3 samples can be made by mixing events. 

We have analysed the three-particle correlations as a function of the kinematical variable 



Q = V 9i2 + 9i3 + ?23 with q^=-{p,-pjf (28) 



Fig P shows -R3, the genuine three-particle correlation function for like-sign triplets, as 
approximated in Eq (0). A strong correlation is observed for small Q- values. There is a 
dip in the curve for Q-values around 1 GeV which is compatible with the depletion of p°'s 
around its mass and gives an indication of the error from using unlike- charged pions in the 
approximation of R3. 
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Figure 8: The Q-dependence of the genuine three-particle correlation function R3, defined in 
the text. 
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